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ABSTRACT 

We announce the discovery of a new early- type, double-lined spectroscopic binary in 
the LH 54 OB association in the LMC. We present a V light curve and radial velocities. 
We investigate the possible configurations of the system, concluding that it probably 
contains the most massive star measured at the present, with a mass of the order of 
100 Mq, while its companion has approximately 50 Mq. 
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1 INTRODUCTION 

Though we can trust that the present development of stellar 
interiors modelling reproduces qualitatively the consecutive 
stages along the evolution of stars, there still remain appre- 
ciable quantitative discrepancies between theoretical predic- 
tions and empirical evidence, especiallly in that concern to 
the higher mass stars. That is not surprising since, on one 
hand, their modelling is intrinsically more difficult due to 
their more complex evolution, and on the other hand, very 
massive stars are rare because of its shorter lifetimes, so that 
there are scarce observational constraints. 

During the last years, we have performed a campaign 
of observations of early Magellanic stars in order to enrich 
the empirical knowledge of these objects, particularly in low 
metallicity regions. 

One of the new binaries discovered for which its light 
curve was completed is a V ~ 13 mag star in the LH 54 OB 
association (Lucke & Hodge 1970). This star (a = 5 h 26 m 24 s , 
5 = — 67°30'13") is bright enough and has sufficiently big 
radial velocity excursions to be at scope of the 2.15-m tele- 
scope at CASLEOQ (San Juan, Argentina) for radial velocity 
measurements. 

Hill et al. (1994) identified this object as LH 54-425 and 
performed UBV photometry, obtaining V = 13.19 ± 0.01, 
B — V = -0.31 ± 0.01 and U - B = -1.02 ± 0.01. Later 
on, Oey (1996a) obtained V = 13.13 ± 0.007, B — V = 
-0.215L0.014 and U-B = -1.010 ± 0.015, listing this star 
as L54S-4. CCD spectra were also obtained by Oey (1996b) 
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with the ARGUS multifiber spectrograph at the CTIO 4-m 
telescope. The spectral type assigned by Oey to this star 
was 04 III (f*). She did not detect double lines at the phase 
of her observation. 

LH 54-425 is the earliest star of the OB association LH 
54, located at the East-side of the superbubble DEM 192 
(Davies, Elliot & Meaburn 1976) (or N51D, Henize 1956) 
that is probably related with its growth (Oey & Smedley 
1998). 



2 OBSERVATIONS 

LH 54 had been included in one of our selected fields for 
searching new eclipsing binaries. It was observed during four 
observing runs between 1998 and 2001. 

LH 54-425 showed no eclipses, but exhibited small peri- 
odic luminosity changes of the order of ~ 0.1 mag, suggest- 
ing an ellipsoidal nature of the variability. For this reason 
we included it as subject for further spectroscopic study as 
soon as its light curve was sufficiently complete to compute 
a reasonably reliable ephemeris. 

2.1 Photometry 

CCD images of the region of LH 54 were obtained at 
CASLEO between 1998 and 2001, with the same instru- 
mentation described in Ostrov et al. (2000). We derived 
a V light curve from aperture photometry performed with 
DAOPHOT (Stetson 1987, 1991), and tied all the measure- 
ments to a unique instrumental system using several stars as 
local standards to reduce errors, as described in Ostrov et al. 
(2000). We performed absolute photometry during a photo- 
metric night (Dec 3, 1998), obtaining V = 13.05 ±0.015, and 
B -V = -0.121 ± 0.011. Table [j] displays our photometric 
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measurements of LH 54-425, corrected to the standard sys- 
tem, together with the internal errors derived from the local 
standards, the FWHM and the airmass for each frame. 



2.2 Spectroscopy 

Spectroscopic observations were obtained by means of a RE- 
OSC spectrograph in its single dispersion mode, attached to 
the 2.15-m telescope at CASLEO. The CCD detector was 
the same Tektronik used for photometry. We used a 600 
mm -1 grating giving a reciprocal dispersion of 1.63 A/pixel 
on the range from 3900 to 5500 A. Our resolution, measured 
from the Cu-Ne-Ar comparison lamp lines, was 2 pixels. 

We obtained spectra near the quadratures in series of 
three 1200 s exposures. The spectrograms were processed 
and extracted using IRAF0 routines. Each series of three 
observations was combined in one spectrogram, before per- 
forming the concomitant measurements. 



Table 2. Radial velocities 
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spectrograms and severely blended in others. The He II ab- 
sorption pairs were measured by means of a double gaussian 
fit using the "deblend" function of the "splot" IRAF task. 
Only the easily separable He II pair s were used for the radial 
velocity determinations. Table 3. J lists the heliocentric ra- 
dial velocities derived for both binary components, together 
with that corresponding to the nebular [O ill] A 5007 emis- 
sion, used to check the stability of the system. 



3 ANALYSIS 

3.1 Ephemeris 

We applied to the photometric data two period search meth- 
ods: one derived from the Lafler & Kinman's (Lafler & Kin- 
man 1965) and the one described in Shwarzenberg-Czerny 
(1997). Because of our short time base and the intrinsically 
small light variations, we can not derive a precise ephemeris. 
We found that the most probable period is P = 2.2475 
days, with aliases each 0.007 days. We used this period and 
Eo = 2452216.597 to plan the spectroscopic observations. 

3.2 Spectroscopic classification 

We classified the spectra following the criteria of Walborn & 
Fitzpatrick (1990). The earliest component shows no lines 
of He I, corresponding thus an 03 spectral type. The other 
component shows He I 4471 hardly visible due our poor sig- 
nal to noise ratio (in fact, it is only detectable in our best 
spectrograms, obtained during dark nights) and more in- 
tense He II lines which is in agreement with an approxi- 
mately 05 type. With regard to the luminosity classes, we 
combined all the spectrograms acquired in each quadrature 
to increment the signal to noise ratio. In the averaged spec- 
trograms, N IV 4058 emission from the primary is visible, 
together with Si IV AA 4089 or 4116 depending on the phase 
(this is alternatively blended with the H<5 absorption of the 
secondary). We are not able to assign a luminosity class to 
the secondary. So, we classify the system as 03 III(f*)+05:. 
Fig. [l] shows the CCD spectrograms of LH 54-425 obtained 
during dark nights, corresponding to both quadratures. 

3.3 Radial velocity measurements 

The He II 4686 line appears clearly double in all our spec- 
trograms. He II 4200 and He II 4542 appear double in some 

2 IRAF software is distributed by NOAO, operated by AURA 
for NSF. 



3.4 Geometry system's exploration 

Since the system shows no true eclipses, it is not possible to 
derive reliable physical parameters just from the light curve 
and radial velocities. For this reason, we assumed an abso- 
lute magnitude and investigate which system configurations 
were compatible with the observations. 

We computed the models with the Wilson-Devinney 
codes (Wilson 1990, Wilson & Devinney 1971). For the 
bolometric albedos and gravity darkening coefficients we as- 
sumed values of A — 1.0 and g = 1.0, respectively, which are 
adequate for radiative envelopes (Rucinski 1969, Lucy 1976). 
For the limb darkening a square root law (Dfaz-Cordoves & 
Gimenez, 1992) was used, taking the corresponding coeffi- 
cients from tables by Dfaz-Cordoves et al. (1995). 

3.4-1 Detached configurations 

We first explored the possible detached configurations, so 
the mode of operation 2 of the Wilson-Devinney program 
was chosen. In order to reduce the degrees of freedom of the 
problem, we fixed the total bolometric magnitude of the sys- 
tem. We considered three cases, with A'hoi ~ —10, —10.3 and 
— 10.6. For the models with Mtoi = —10 we assigned temper- 
atures of 52500 K and 44500 K, corresponding to spectral 
types 03 V and 05 V according with Schmidt-Kaler (1982), 
while for those models with Mboi = —10.6 we adopted 
the temperature scale of Chlebowski & Garmany (1991) for 
03III and 05III stars, 46500 K and 42300, respectively. For 
the models with intermediate bolometric magnitude we used 
the mean of the two temperature scales. With this proce- 
dure, the extreme cases in distance moduli and temperature 
scales were considered. 

For each adopted orbital inclination, temperatures and 
potential fii of the primary, the secondary's potential O2 
was chosen according with the assumed total bolometric 
magnitude of the system. 

Fig. [2] displays the results of the experiments with de- 
tached configurations. Points of bigger size represent smaller 
O-C light curve residuals. The lower limit corresponds to the 
primary being 2.5 magnitudes brighter than the secondary 
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Table 1. V light photometry of LH 54-425 
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component, which can be rejected according to the relative 
intensity of the corresponding spectral lines. The upper limit 
corresponds to semi-detached configurations, with the sec- 
ondary star filling its Roche-lobe. Models above the solid 
curve have the secondary brighter than the primary, there- 
fore can be also disregarded. 

The upper panel, that corresponds to highest tempera- 
tures and lower luminosity, shows that the best fitting mod- 
els are near the semi-detached configuration, with the sec- 
ondary filling its Roche-lobe. These models are not realistic 
from an evolutionary perspective, since they have secondary 
components of bigger size than the primary. If the less mas- 
sive star is the more evolved one, we have to conclude that 



mass inversion has occurred, and the system is now semi- 
detached. 



In the middle panel we considered models with inter- 
mediate luminosities and temperatures. The best solutions 
are aligned at the upper limit of the models considered, in- 
dicating a semi-detached configuration. There are also good 
models with i ~ 58deg and fii ~ 3. It means that the pri- 
mary star is almost filling its Roche-lobe. However, it is not 
very plausible since the magnitude difference between both 
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stars is near 2.5 mag, which does not agree with the spec- 
troscopic evidence []. 

The residuals from models with lowest temperature 
scale and highest total luminosity (which allows the biggest 
sizes of the stars) are presented in the bottom panel. In this 
case, besides of the semi-detached configurations, there ex- 
ist plausible detached solutions, with i ~ 47 deg. The lowest 
value considered for the potential fii (2.85) implies that the 
primary component is almost filling its Roche-lobe, although 
(as mentioned above) in these models the luminosity differ- 
ence between both stars is excessive. 



3.4.2 Semi-detached configurations 

Considering that some results of the experiments suggested 
that the system is semi-detached, with the cooler and less 
massive star filling its Roche-lobe, we explored these con- 
figurations, using therefore the mode of operation 5 in the 
Wilson-Devinney programs. Figure ^| shows the residuals of 
these models, corresponding to the high and low tempera- 
ture scales above mentioned, respectively. As the margins of 
plausible models are wide, it is necessary to consider limits 
to the total luminosity to restrict the amount of solutions. 
The upper and bottom curves drawn on figure ^ correspond 
to models with Mboi — — 10 and Mbol = —10.6, respectively, 
while models above the horizontal line have the secondary 
more luminous than the primary component. Configurations 
with Q\ smaller than those plotted in the graphics are in 
overcontact, and can be also rejected according to the as- 
sumed limits in the total luminosity of the system. 



3 It does not allows to disregard completely such solutions. 
For example, Niemela & Bassino (1994) called "primary" at the 
Roche-lobe filling component of HV2543, based on the intensity 
of their spectral lines, while from the W-D modelling that star is 
approximately 0.9 V magnitudes fainter than its companion. 



4 RESULTS 

In spite of the wide range of models that can fit well the 
observations, we can state with reasonable confidence that 
the inclination is roughly 45 ~ 50° , since bigger values would 
cause eclipses. This implies that the masses are ~ 100 and 
~ 5OM0 for the 03 and the 05 component, respectively. 
Table H summarises the parameters that best fit suitable 
detached and semi-detached solutions. It must be noted that 
the errors quoted in Table ^ correspond to fixed values of 
i, fii and £^2, i.e., the errors of the masses reflect only the 
errors in the radial velocity measurements. 

The modelled light and radial velocity curves are dis- 
played in Fig. [I and Fig. ^, together with the observations 
and O-C residuals. The models correspond to the best fit- 
ting detached configuration, although there are not visible 
difference with the best fitting semi-detached solution. 

In the more modest case, the inclination is i ~ 58° and 
the corresponding masses are ~ 64 and ~ 31Mq, although 
this model have a primary 2.3 magnitudes brighter than the 
secondary component, that is not supported by the spectro- 
scopic data. (The other extreme case, in that both stars have 
similar brightness, corresponds to i ~ 39° and gives masses 
of ~ 160 and ~ 80Mq). All possible solutions indicate that 
this system contains the more massive star measured at the 
present. 

In a previous work (Ostrov 2001) we found from the 
analysis of another early Magellanic binary with 03If* and 
O6V components, masses of 41 and 27 Mq. That was an 
overcontact system, with a still shorter period, of only 1.4 
days. These results suggests that extremely early stars have 
a wide range of masses, at least when they belong to very 
close binaries, where mass exchange and mass loss play a 
fundamental role in the evolution. 

Indeed, observations of these systems with higher signal 
to noise and dispersion than that reachable at CASLEO 
would provide valuable data to investigate the similitudes 
and differences between these objects. 
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Figure 3. Residuals of the semi-detached models. Points of differ- 
ent sizes represent the sums of squared residuals, as in Fig. ^. Top 
and bottom panels correspond to different assumed temperature 
scales, as quoted in the text. 
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Figure 2. Residuals of the detached models. Points of different 
sizes represent the sums of squared residuals, the bigger symbols 
corresponding to the smaller amounts. The limits between dif- 
ferent point sizes are 1.65, 3.30, 4.95 and 8.24 times the sum of 
squared residuals of the best fitting model. The different panels 
correspond to different assumed effective temperatures and total 
luminosity of the system, as quoted in the text. 
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Figure 4. The observed and modelled V light curve for LH 54- 
425 and their concomitant O-C residuals. 
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Figure 5. Observed and modelled radial velocity curve for LH 
54-425 and their O-C residuals. Hollow circles correspond to the 
primary component and filled ones stand for the secondary. 

Table 3. Best fitting plausible models 





Adopted Values 








46500 K 


T 2 




42300 K 


M boi 




-10.6 




Fitted Values 




V 7 




295.0 ±3 km s _1 


q (Mi/Mi) 




0.48 ± 0.01 


configuration 


detached 


semi-detached 


i (adopted) 


47° 


45.5° 


Hi (adopted) 


3.1 


3.15 


n 2 


2.86 a 


2.84 b 


a 


38.1±O.4R 


39.3 ± 0.4 R 


Eres 2 


0.0033222 


0.0033563 




Star Dimensions 




Mi 


100 ±3 M 


108 ±3 M 


Ri 


15.1 R© 


15.1 R 


M hol i 


-10.2 


-10.2 


loggings] 


4.1 ± 0.1 


4.1 ±0.1 


M 2 


48.4 ±3 M 


52.2 ±3 M 


R2 


11.9 R Q 


12.4 R 


M bol 2 


-9.3 


-9.3 


l°g92[cgs] 


4.0 ±0.1 


4.0 ±0.1 



a Chosen according the total luminosity of the system. 

k ^crit 



